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Abstract 



We present a mean field solution of the antiferromagnetic Heisenberg chain 
with nearest (Ji) and next to nearest neighbor (J2) interactions. This so- 
lution provides a way to estimate the effects of frustration. We calculate 
the temperature-dependent spin- wave velocity, v s (T) and discuss the possi- 
bility to determine the magnitude of frustration J2/J1 present in quasi ID 
compounds from measurements of v s (T). We compute the thermodynamic 
susceptibility at finite temperatures and compare it with the observed sus- 
ceptibility of the spin-Peierls compound CuGeC>3. We also use the method 
to study the two-magnon Raman continuum observed in CuGeC>3 above the 
spin-Peierls transition. 
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I. INTRODUCTION 



The discovery of the spin-Peierls (SP) transition in the inorganic material CuGe03 [JTJ] 
has led to an intense investigation of the magnetic properties of this system. It is becoming 
increasingly clear that this material may not be a prototype SP system such as the quasi one- 
dimensional organic material TTFCuBDT ||. For example, the temperature dependence of 
the SP gap || is unlike what is expected for conventional SP systems 0. A recent inelastic 
neutron scattering experiment reports observing a spin gap at temperatures above the SP 
transition temperature T$p. So far, there has been no evidence for the presence of a phonon 
soft mode ||. X-ray scattering studies of the incommensurate phase of CuGe0 3 show the 
existence of a soliton lattice with the width of the soliton being much larger than that 
predicted by calculations @. These results, taken together, call for a better understanding 
of the homogeneous state of CuGe03 above Tsp which should eventually shed light on the 
nature of the SP transition itself. 

The basic structure of CuGe0 3 consists of edge-sharing Cu0 6 octahedra forming Cu0 4 
chains along the crystallographic c axis. The dimerization of the S = | Cu ions below 14 
K has been determined by neutron- diffract ion measurements [0]. From inelastic neutron 
scattering (INS) measurements, Nishi et al. || estimated the intrachain and interchain 
exchange parameters J c ~ 120 K, Jf, « 0.1J C , suggesting that interchain effects may not 
be negligible in this compound. On the other hand, INS measurements above T$p 
clearly show the two-spinon continuum which is characteristic of a one dimensional 5* - 
antiferromagnetic chain []10[. The susceptibility of CuGeOs above T$p M shows a broad 
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maximum as expected for Heisenberg chains. From the temperature at which this maximum 
is observed (56 K), one can estimate the value of J c using the results of Bonner and Fisher 
||1 1|| , yielding J c = 88 K. This discrepancy between the value of J c estimated from static 
susceptibility measurements and INS led to the proposal [0, [|l3j that a minimal model to 



describe the magnetic properties of this system above T$p is the so-called "Ji - J2" model. 
This model describes an antiferromagnetic Heisenberg chain with nearest neighbor (n.n.) 
and next-to-nearest neighbor (n.n.n.) exchange interactions. The model Hamiltonian is 
written as 

H = J^2(Si ■ Si+i + aSi ■ S i+2 ) , (1) 

i 

where J = J c is the intrachain superexchange between neighboring Cu ions along the c 
direction. The second term in the Hamiltonian ([I]) is the exchange interaction between next 
to nearest neighbor Cu ions. In CuGeOs, the n.n.n. superexchange path is through Cu - O 
- O - Cu and is identical to that in the cuprate superconductors. A detailed analysis of the 
structure of CuGeOs and its relation to the magnetic interaction can be found in |14],15 



The model Hamiltonian (JJ) has been studied by several authors [^6|. Though the Hamil- 



tonian is not exactly solvable for all values of a, the phase diagram is well understood qual- 
itatively. For < a < a cr , the ground state remains gapless (as is the case when a = 0). 
The effect of a is to renormalize the spin wave velocity in this regime. The value of a cr has 



been estimated to be 0.2411 [17]. When a > a cr , the spectrum becomes gapped and for 



a — 2, the Hamiltonian (|XJ) is exactly solvable [I8|. As mentioned earlier, there has been a 



renewed interest in this model in the context of CuGeOs. R-iera and Dobry [|T^] as well as 



Castilla and co-workers |13| computed the thermodynamic susceptibility of the Hamiltonian 
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(H) numerically and compared it with the experimental values. Both groups found that the 
presence of a non- vanishing a is needed to provide a consistent description of both the INS 
and susceptibility results. 

If indeed the J\ - J 2 model is an appropriate starting point to embark on a study of 
CuGeC>3, it would be desirable to have some way of calculating physical quantities, especially 
considering the wealth of experimental data now available on CuGeC>3. This is the primary 
objective of this study where we present some results obtained from a solitonic mean field 
theory of the J\ - J2 model. This method, which is meaningful only for a < a cr , provides 
a simple self-consistent way of evaluating the effects of frustration. We use this solution to 
obtain the spinon dispersion relation as a function of a. We also calculate the ground state 
energy and bulk susceptibility at finite temperatures as a function of a. The results for 
the bulk susceptibility are compared with the experimentally observed values in CuGeO"3. 
We then examine the Raman continuum seen experimentally above T$p and show how this 
arises as a natural consequence of competing magnetic interactions. We compute the Raman 
intensity using the mean field solution and compare with experimental results. The paper 
is organized as follows. In section II, we present the mean field solution of the Hamiltonian 
(H) and use it to compute the static susceptibility. In section III, we compute the two- 
magnon Raman scattering intensity arising from competing magnetic interactions. Section 
IV contains a brief summary of our results. 



II. A MEAN FIELD SOLUTION OF THE J x - J 2 MODEL 

In this section, we propose a mean field solution of the Hamiltonian ([I]). This solution is 
based on a mapping introduced by Gomez-Santos |1| between the spin-^ Heisenberg chain 



with n.n. interactions and a Hamiltonian describing the dynamics of antiferromagnetic 
domain walls. It has also been used by Weng and collaborators |2(| for the one dimensional 



t - J model. We show below how the mapping can be generalized to the spin Hamiltonian 
(H) with n.n.n. interactions as well. In this mapping, the local degrees of freedom are given 
by the nature of the bond (ferromagnetic or antiferromagnetic) between two interacting 
spins. We work using periodic boundary conditions and choose the following convention: 
Neel ordering is characterized by an "up " spin at the first site. Since, by definition, the Neel 
ordered state does not have a "kink" (henceforth, we shall use the words "kink", "soliton" 
and "domain wall" interchangeably), this state is the vacuum state |0) of the solitons, which 
we write symbolically as 

|0> = I TITITITITITI • • •) • 

A state with a kink between sites i and % + 1 is defined to be a one-soliton state d\\0). For 
example, d\\0) defines the spin configuration 

4io> ee 1 mmutm-) • 

With these definitions, it is easy to verify that the n.n. term in (|l|) is mapped into 

- £ {(4-Wi+i + 4-i4+i + h.c.)(i - d\d t ) + (djdi - i)| . 
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Here, we have avoided the sign-problem by assuming the solitons to be hard-core bosons. 
On performing the Jordan- Wigner transformation 

di -> exp (z'7T dlndm) di , 

one sees that the n.n. term written above preserves its form. By looking at the action of 
the n.n.n. term in ([!]) on a pair of interacting spins (or equivalently, a given bond), one can 
write down the n.n.n. term in terms of the soliton operators in the same manner as above. 
Doing this and performing a Jordan- Wigner transformation, we obtain the J\ - J 2 model in 
terms of fermionic soliton operators as 

H = \ E {(4+i + - 4^)(4-i - *-0 + (4di - i)} 

+ -7T E(4+2 + di +2 ){d\ +1 di + 4<W(4-i - ( 2 ) 
^ i 

+ ^^(24rf,-l)(24 +1 rf m -l) . 

i 

The above mapping may also be verified by using the definition of the original spin operators 
in terms of the solitonic operators, viz., 

\{4-i- d i _ l ){4 + d i ){i-2s*) 

\{4-i-d i - l ){d\ + d i ){l+2S* i ) 

U-y+ l e^^Y.d)d 3 ) . 

Substituting the above expressions for the spin operators in (|l|), we recover (0). It should 
be noted that the maximum number of fermion operators occurring in (Q) is four, both for 
the n.n. term and for the n.n.n. term. Longer ranged interactions such as S« • Sj + 3 would, 
on the other hand lead to terms containing products of six or more fermion operators. 

The Hamiltonian (|2|) is solved by treating the quartic terms in mean field theory. We de- 
fine the following averages that are determined self consistently: n = (d\di), Aj = (<ij_i^i+i) 
and A 2 = (4-i4+i)- m terms of these averages, the mean field Hamiltonian is given by 

Hmf = Y,dUi d[l - 2(Ax + A 2 )] + aJ(2n - 1)} 
i z 

+ E(4-i^+i + 4-A+i + h.c.) {^(1 - n) + aJ(Ai + A 2 )} 
+ E , 

where N is the number of spins in the chain and Eq is defined by 
E = jVJ(Ai + A 2 ) - - NJ(1 - n)(A 1 + A 2 ) - NaJ^ + A 2 ) 2 - NaJ(n 2 - J) . 



s+ = 
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The mean field Hamiltonian can be simplified further by defining the following two quantities 

Jb = ^[1 - 2(Ai + A 2 )] + aJ(2n - 1) 
J A = ^(l-n) + aJ(A l + A 2 ) . 

Using these definitions and Fourier transforming, the mean field J\ - J 2 model is written as 
H MF = ^2{Jb + 2J A cos2k)d\d k + J A wa.2k{d\S_ k - d_ k d k ) + E , (3) 

k k 

where the lattice constant has been set to be unity. The mean field Hamiltonian (^) can 
now be solved by introducing the Bogoliubov transformation 

d k = u k a k - iv k a ] _ k , 

where 

Vk= V5 ( E~k )2 S&nk ' 



with e k = Jb + 2J A cos2fc, A^ = 2J A sin2fc and E k = y£ k + A£. It is then easy to see that 
the Hamiltonian (|3]) reduces to 

H M f = E k®k a k + \ J2( e k - E k ) + E . 

k k 

On evaluating the mean field quantities and substituting them in the definitions for J A and 
Jb we get two mean field equations that have to be solved for self consistency: 

J 1 ^ J B cos 2k + 2 J A a J ^ 2 J A cos 2k + J B _ , 
Jb = 2 [1 + aF E ^ (1 - 2«t)l " 77 E ^ (1 - (4) 

27, = ^ [1 + 1 S 2J ^ + A (1 - 2n t) ] - % £ + _ ^ 

where n*. is the usual Fermi distribution function. When a — 0, the above equations reduce 
to the mean field equations written down in [B5 and EO]. The solution in that case was 



obtained as Jb = 2J A . By inspection, we see that the solution of the mean field equations 
is still given by Jb = 2J A = J, where now, J is modified by a. The dispersion relation is as 
before, E k = 2J\ cosfc|. On substituting the solution in either of the above two equations, 
we determine J by the equation 



1 — 2a 

1 H — — I cos k\ tanh(/3J| cos k\) 



N k 
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At zero temperature, the above expression gives us the dispersion relation, 



E k = J[l + ^ — ^1] | CO sA;| . (5) 



7T 



When a = 0, @ describes the dispersion relation for spinons obtained by Faddeev and 
Takhtajan As pointed out in ||20|| , the spinon velocity obtained from the mean field 



theory of antiferromagnetic domain walls, namely (1 + 2/tt) J ~ 1.64 J is quite close to 
the exact value (ir/2J « 1.57 J). When a is nonzero, our result shows that the spinon 
velocity is reduced because of frustration arising from the n.n.n. interaction. This is to 
be expected from physical grounds. On comparing the spinon velocity in (^) with the 
spinon velocity for a = 0, we see that the effect of the n.n.n. term is to reduce the spinon 
velocity as v s (a) = v s (0)[l — 4a/ (tc + 2)]. The ground state energy at T = is given by 
E = E (a = 0) + 2aNJ/ir 2 , where E (a = 0)/NJ = - (1/vr + 1/tt 2 ) « -0.420 is the 
ground state energy of the n.n. Heisenberg chain within this approach, which compares well 
with the exact result 1/4 - In 2 « -0.443 |22|]. 



The mean field solution described above provides a simple way to calculate physical 
quantities in the J\ - J2 model. In particular, one can see how the presence of the n.n.n. 
term in the Hamiltonian ([!]) alters the susceptibility at finite temperatures. Before we 
proceed to show this, we discuss the limitations of the mean field theory. First, we point 
out that though the mean field solution is valid for any value of a, the theory itself is 
meaningful only in the gapless regime, i.e., when a < a cr where the spinons are deconfined. 
This mean field theory is not suited for studying the transition between the gapless and 
gapped phases of the J\ - J 2 model. We reemphasize that our objective is not to undertake 
such a study (indeed, far more powerful techniques are available and have been used to 
study this problem) but to provide a simple self consistent prescription to calculate physical 
quantities in the gapless phase of the J\ - J2 model. Next, we point out that the mean field 
solution is plagued by its reluctance to move away from the universality class of the XY 
model. This can be seen by writing down the (mean field) ground state wave function in 
terms of the soliton operators d k as \^) = l\ k (u k + Vkd k d_ k ) |0). Whereas in the XY model, 
the Hamiltonian can be solved exactly to determine u k and v k , the mean field solution of 
the Heisenberg as well as the J\ - J 2 model give the same u k and v k . In this sense, this 
theory has the same difficulties as the mean field theory of Bulaevskii ]23[] though the nature 
of the quasiparticles in the two theories are entirely different. Finally, let us consider the 
reduction in the spinon velocity as a function of the n.n.n. coupling a. From flg), we see that 
the velocity v s (a) ~ v s (0)(l — 0.8a). This result can be checked numerically. In a recent 
numerical study, Fledderjohann and Gros |24| found that the spinon velocity in the J\ - J2 
model can be fit to the relation v s (a) = v s (0)(l — 1.12a). Thus we see that the results from 
mean field theory underestimate the effects of frustration. 

We are now ready to consider the static susceptibility at finite temperatures. In particu- 
lar, we are interested in seeing how the presence of a n.n.n. interaction alters the susceptibil- 
ity of an antiferromagnetic Heisenberg chain. Now, it is well known that the susceptibility 
of the Heisenberg chain is constant at zero temperature and increases to a broad maximum 



at T = 0.641J [|TTJ. The increase in susceptibility is due to the gapless spinon excitations. 
As we shall see, the n.n.n. term causes a suppression in the maximum value of this sus- 
ceptibility. The reason for this is that the spinon velocity now has a different temperature 
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dependence with a nonzero a. To illustrate this, we show in Fig. p], the spinon velocity as 
a function of temperature. The long-dashed line shows the results for the case a = 0. The 
solid line shows the spinon velocity as a function of temperature for a = 0.24. For the latter 
case, we have, for purposes of illustration, chosen J such that the velocity at T = 0, given 
by (H) is the same as that for a = 0. It should be noted that while the low temperature 
velocities for different a are identical, the results begin to differ as temperature increases. 
In particular, they are perceptibly different at T « 0.6J. Therefore, we should expect a 
difference in the finite temperature susceptibilities as a function of a at these temperatures. 

We have calculated the static susceptibility at finite temperatures from the mean field 
solution. At this point, it is crucial to realize that the results from the mean field theory 
violate spin-rotational invariance of the underlying Hamiltonian, and that the physically 
relevant susceptibility would be given by the average 

_ 1, 

X g {Xxx + Xyy "I" Xzz ) i 

where Xxx, Xyyi Xzz are the susceptibilities for the applied magnetic fields in x, y and z 
directions respectively. Next we note that our mean field theory becomes exact for the XY 
model for which the susceptibilities x^I = 1/vr ~ 0.318 and Xxx = Xy V Y ~ 0.075 « x^I '/4 



25| , in dimensionless units and at zero temperature. We therefore have for the XY model, 



the relation 

X « *£72 • (6) 

We have applied this relation to our mean-field solution of the Heisenberg model. 

A straightforward calculation shows that the susceptibility at zero temperature is given 

by 

X 1 
g^ 2 B N 2nv s ' 

where g is the gyromagnetic ratio, fiB is the electron Bohr magneton and v s , as before, is the 
spinon velocity. The factor 2 in the denominator comes from the average over the directions 
of the applied magnetic field as given by Eq. (||). 

As a quick check of our results, let us consider the case a = and compare with the results 
of Bonner and Fisher ||11|| . At T = 0, from (|5|) and the above expression for susceptibility, 
we get x{T = 0) J/(g 2 fi 2 B N) « 0.106, which compares very well with the Bethe-Ansatz 
results 1/tt 2 ^ 0.101 |TT]]. At finite temperatures, we solved the mean field equations (f|) 
numerically and used it to determine the susceptibility. We found that the temperature at 
which the susceptibility attains its maximum value is given by fcsT max /J = 0.61, while the 



corresponding value quoted in [J Ij is 0.641. These comparisons encourage us to proceed to 
the case a^O. 

Our results for the case are summarized in Fig. ^] where they are compared with 

experimentally obtained values from a polycrystalline sample of CuGe03 [p6| |. We have 



chosen g = 2.0 in all our calculations. The long-dashed line shows the results for J = 88 K, 
a = that are clearly in disagreement with experiment. The choice of J = 88 K stems 
from the observed temperature (56 K) at which the susceptibility attains its maximal value 



and its relation to J from the results of [11 1 or equivalently, our results. The dotted line 



shows the results we obtained from choosing J = 115 K, a = 0. This choice of J yields the 
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observed value of spinon velocity from INS [§] without including the effects of frustration. 
It is clear that there still remains a discrepancy between the observed and calculated values 
at intermediate temperatures. Our results for the susceptibility with a non vanishing a are 
given by the solid line. Here, we have chosen J = 142 K and a = 0.24. This choice of 
parameters when substituted in (^]) gives the value of the spinon velocity observed in INS. 
(Note that the corresponding value of J determined numerically in |13| is J = 150 K.) At 



low temperatures, results with and without a are identical (see Fig. [I| and our discussion of 
those results). However we see that at as temperature increases, a non vanishing a does tend 
to suppress the susceptibility. We see that there is still a discrepancy between theory and 
experiment which is of the order of 10%. We believe this discrepancy is because the mean 
field theory underestimates the effects of frustration. We also considered the possibility that 
the discrepancy between theory and experiment may well be due to the effect of interchain 
coupling. As mentioned earlier, the interchain superexchange Jb ~ 0.1J C and this could play 
an important role. A qualitative argument against this possibility is that at temperatures 
where the discrepancy is maximal (T > 0.5Tsp), the magnetic correlation lengths should be 
small. To substantiate this argument, we did a simple calculation using renormalized spin 
wave theory with the appropriate values of Jb and J c . Were it to order antiferromagnetically, 
the observed values of the superexchange interactions suggest that the "Neel temperature" 
of CuGeOs would be ~ 10 K which is lower than T$p. At temperatures above T$p, the effect 
of (short range) antiferromagnetic ordering can be easily estimated by spin wave theory. A 



standard calculation using the "decoupling approximation" |27j gives us the susceptibility 
X in terms of the integral equation 

1 f w / >7r , , 1 1 

— / dz dy- 



7r 2 Jo Jo " 2(7 -Ik) -X 1 4k B T 

where 7 = J c + Jb and 7^ = J c cos;z + Jbcosy. We have ignored J a ~ — 0.01J C and set g\ip 
to be unity for convenience. One integration in the l.h.s. of the above expression can be 
done trivially and we obtain 

1 7T / J 

dz- 



J(X + cos zf -0.01 2 \k B T 



where x = {2JX) 1 — 1-1 and J c = J = 10 Jb- Solving the above equation numerically 
with J = 115 K, we obtain the susceptibility as shown in Fig. |[ It should be remembered 
that the theoretical results are valid only at temperatures well above T$p where the effects 
of phonons as well as effects due to competing antiferromagnetic and SP instabilities can 
be ignored. The results in Fig. ||] and Fig. [| show that a one dimensional model with 
frustration is indeed a better starting point to study CuGe0 3 rather than an anisotropic 
two dimensional model. This is in agreement with a number of experimental observations 
such as (i) a non- zero nuclear spin-lattice relaxation rate seen in a Cu NQR study p8|] , (ii) 
the two-spinon continuum observed in INS above T$p 0,0, and (iii) a broad continuum of 



magnetic excitations seen in two-magnon Raman scattering above T$p P5 |, which will be 
discussed in the next section. All these experiments show characteristic features of a one 
dimensional Heisenberg antiferromagnetic chain. 

Since our mean field solution underestimates the effects of frustration, we are unable to 
determine the exact ratio of J2/J1 in CuGeOs. As we mentioned in the introduction, the 
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value of Jij J\ was estimated to be 0.36 Jl2| and 0.24 |13j by two independent studies. Since 
these two values lead to two different fixed points (gapped and ungapped spin excitation 
spectra respectively), careful experimental investigation of spin dynamics above T$p may 
provide further clues. Recently, Kuroe et al. have claimed that J 2 / ' J\ = 0.35 ± 0.05 
leads to a good fit of the magnetic specific heat as determined by quasi elastic Raman 
scattering. It should be noted, however, that this choice does not fit the observed magnetic 
susceptibility. CuGeC>3 is 0.35 ± 0.05. 

Before we conclude this section, we discuss how the temperature dependence of the 
spinon velocity may be used to determine the presence of competing magnetic interactions 
in quasi ID materials. As shown in Fig. [1], a non zero value of J2/J1 leads to a difference 
in the temperature dependence of v s (T). This feature may be detected experimentally. 
As a possibility, we suggest INS experiments in the presence of a magnetic field. In the 
presence of a magnetic field, it is known that the spinon dispersion goes to zero at k 
given by cos A; = —2H/v s (T) f3"I|| . This is a feature of low dimensional spin excitations 
and is seen in the XY model as well. This behavior should be contrasted with spin wave 
(magnon) excitations, where no such node appears in the magnon dispersion as a function 
of the magnetic field. An examination of the temperature dependence of this node with 
moderately high fields H/J ~ 0.2 should indicate the presence of frustration as shown in 
Fig. 0- 



III. TWO-MAGNON RAMAN CONTINUUM IN CUGEO3 

In this section, we consider the Raman spectra in the homogeneous phase of CuGeC>3 
observed in two-magnon scattering experiments ||29|| . A partial account of the results in this 
section has been published elsewhere [|3"2]| . For completeness, we first give a brief summary 



of the pertinent experimental results. 

Scattering of light from magnetic (as opposed to phonon) excitations in CuGeO"3 is 
strongest in the (zz) geometry viz., when both incident and scattered light are polarized 
along the crystallographic c direction. This is consistent with the fact that exchange in- 
teractions are strongest along this direction. The Raman spectra above T$p are rather 
featureless and the spectral weight is distributed over a wide range of energies (100-300 
cm -1 ). A broad maximum is observed around 250 cm -1 . This should be contrasted with 
the Raman spectra observed in antiferromagnetically ordered compounds such as the ru- 
tile or perovskite halides. In these compounds, one typically observes a well defined peak 
at characteristic magnon energies which is broadened by magnon- magnon interactions |[33|| . 
This comparison immediately suggests light scattering in CuGeO"3 from a continuum of 
excitations rather than single-particle excitations (magnons). 

The magnetic Raman scattering in CuGeOs can be studied using the scheme of Fleury 
and Loudon |33| . The basic idea is that in an exchange-coupled system, an incident photon 



can create particle-hole excitations accompanied by a pair of spin flips that propagate as 
magnons. In this scheme the Raman operator Hr has essentially the same form as the 
Heisenberg exchange interaction, i.e. Hr ~ Yli Si • Sj+i. In one dimensional chains with n.n. 
exchange interaction, we are faced with an intriguing situation where the Raman operator 
for scattering in the (zz) geometry is proportional to the Hamiltonian itself. Thus, there 
should be no Raman scattering from a n.n. Heisenberg spin chain, a fact which has not been 
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appreciated previously [p9| |. However, there are two ways in which Raman scattering can 
occur: (i) interchain couplings and (ii) exchange interactions that are extended in space, i.e. 
n.n.n. interactions. Let us first consider the former possibility. The role of the interchain 
couplings is to establish short range antiferromagnetic order. In such a case, incident light 
is scattered off paramagnon excitations. However as we mentioned earlier, the featureless 
spectra observed in CuGeC>3 seem to indicate scattering from a continuum of excitations. 
This is basically a one-dimensional feature. The second possibility is therefore a more viable 
explanation of the experimental results. In view of the results discussed in the previous 
section, it is also a more natural choice. Thus, we work with a Hamiltonian which has both 
n.n. and n.n.n. interactions. This implies that the Raman operator will also have both n.n. 
and n.n.n. terms. For convenience, we subtract from the Raman operator, a part which 
commutes with the basic Hamiltonian of CuGeC>3 (|I|), we can write down a Raman operator 
of the form 

Hr = A ^2 Si • S,; +2 . 

i 

The Raman intensity is then evaluated by considering the appropriate spectral function. 

We now use the mean field theory described in the previous section to compute the 
Raman intensity. The Raman operator R can be written in terms of the soliton operators and 
transformed into quasiparticle or spinon operators just as we did in the previous section. We 
find that there are three basic physical processes caused by R: (i) spinon-spinon scattering, 
(ii) two-spinon creation/annihilation and (iii) four-spinon creation/annihilation. Of these, 
(ii) corresponds to one-magnon excitations. On using the mean field solution, one can show 
that the contribution to light scattering from these excitations vanishes, which is just the 
statement that one-magnon excitations with k = have zero spectral weight. Two-magnon 
scattering is given by (iii). We find that the matrix element for two-magnon (four-spinon) 
creation (Stokes component of the scattering intensity) is given by 



S(k + k' + q + q') a\ a\, at at, 

kk'qq' 



[±exp{-i(2k + k' - q')}{{u k u k > + iv k u k >)(v q v q > + iv q u q >) 

-(u k v k > + iv k v k >)(u q v q > + iu q u q >)} 
- exp{-i(q + q')} u k v k >u q v q > ] , 

where the summation is over the range of spinon momenta and the quantities u k and v k are 
as defined in the previous section. The matrix element can be understood as arising from a 
decomposition of two magnon excitations into four spinon excitations with total momentum 
zero. Given this matrix element, we can calculate the Raman intensity. The mean field 
solution thus provides a way of incorporating matrix element effects arising from the Raman 
operator instead of postulating ad hoc filtering functions that mimic matrix element effects. 
As noted elsewhere ||32|| , such effects could be very important and ignoring them often leads 
to erroneous conclusions. 

The results for the Raman intensity at T = 20 K with J = 142 K and a = 0.24 are shown 
in Fig. |4] (solid line) where they are compared with experimental results (squares) obtained 
from a single crystal sample. Phonon lines at 184 cm" 1 and 330 cm -1 have been subtracted 
from the experimental data. Since we do not know the value of A in our definition of the 
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Raman operator R, we have normalized the maximal theoretical value of the intensity to the 
maximal experimental value. A uniform background of 50 counts has also been subtracted 
from the experimental data. We see that there is a reasonable agreement between theory 
and experiment as far as the observed continuum is concerned. The shoulder observed 
experimentally at ~ 390 cm^ 1 is not understood and our calculations show no indications 
of this feature. We also see that the maximum of the experimental data is shifted slightly 
toward larger energies than the theoretical results. This could either be due to interchain 
couplings and phonon effects that have been neglected or a mismatch in the parameters we 
have chosen. Our calculations are in good agreement with results from exact diagonalization 



as well as those of Singh and collaborators |35] who have studied the same problem 
recently using a finite temperature Lanczos method. We have also calculated the Raman 
intensities at higher temperatures. We find that the continuum broadens out as temperature 
is increased and there is no accumulation of spectral weight at lower energies as seen in 
experiment ||29|| . 

From the results shown in Fig. |4], we conclude that the Hamiltonian (|]) provides a good 
description of the observed Raman continuum owing solely to the presence of the n.n.n. 
interaction. These results as well as the reasoning behind them also lead us to conclude that 
there would be no inelastic Raman intensity in one dimensional chains such as KCuF 3 where 
n.n.n. interactions are negligible. In such materials, while INS can observe the one-magnon 
continuum, the two-magnon continuum cannot be observed in Raman scattering. These 



conclusions are indeed borne out by experiments [36 



IV. SUMMARY 

To summarize the results of this study, we have proposed a mean field solution to the 
antiferromagnetic Heisenberg chain with n.n. and n.n.n. interactions. This was motivated 
by recent interest evinced in the spin-Peierls compound CuGeOs where it is suspected that 
n.n.n. interactions play an important role. Our mean field solution provides a self consistent 
way to include the effects of frustration and calculate physical quantities of interest as a 
function of temperature. Though we have limited our attention to experimental results in 
CuGeOs, the results in this paper can be used to study the effects of frustration in generic 
quasi ID systems. 

As a first application, we have evaluated the temperature dependent spin-wave velocity, 
v s (T) and found that (i) v s (T) decreases rapidly for a = with increasing T, a being the ratio 
between the n.n.n. and n.n. interactions and (ii) this decrease is much less pronounced for 
a 7^ 0. We propose that this fact could be exploited to determine the frustration parameter 
a directly from measurements of v s (T) e.g., by inelastic neutron scattering studies in other 
quasi ID systems. 

We then used the mean-field theory to evaluate the static susceptibility x- A compar- 
ison with experimental results shows that the inclusion of frustration does reproduce the 
suppression of susceptibility at finite temperatures seen experimentally. We have been able 
to identify the physical reason for the suppression of \ as a function of a with the differing 
behavior of the respective v s (T). We have also presented the results of a spin wave cal- 
culation for an anisotropic two dimensional Heisenberg model which do not, on the other 
hand, compare well with the experimental results for the susceptibility. This suggests that 
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a one dimensional model with frustration is a more appropriate starting point to model the 
homogeneous phase of CuGe0 3 . 

We also used the mean field solution to calculate the two-magnon Raman intensity from 
a one dimensional antiferromagnetic chain and showed that frustration is necessary for ob- 
taining a non-vanishing spectral weight. Considering the fact that n.n.n. interactions are 
needed to produce a good fit to the static susceptibility of CuGeOs, using the same idea to 
explain the Raman continuum provides a consistent picture as far as modeling goes. Quite 
apart from the relevance to CuGeO"3, this method of computing Raman intensity can be 
used to study experimental results now available in Sr 2 Cu0 3 and SrCu0 2 |J7J]. Clearly, 
a complete understanding of CuGeOs calls for the inclusion of interchain and phonon ef- 
fects. The mean field solution presented in this paper provides a way to incorporate these 
effects systematically. Work is in progress on these issues and results will be reported in a 
forthcoming publication. 
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FIGURES 



FIG. 1. Spinon velocity as a function of temperature T, obtained by the solitonic mean-field 
theory. The long-dashed line shows results for J = 115 K and a = 0. The solid line is obtained with 
J = 142 K and a = 0.24. The values of J are chosen such that the velocities at zero temperature 
are the same in both cases. 

FIG. 2. Experimental susceptibility of CuGe03 (filled squares, |26fl ) as a function of tempera- 
ture compared with results from the mean field theory of the J\ - J2 model. The long-dashed line 
shows results with J = 88 K and a = 0. The dotted line shows results with J = 115 K and a = 0. 
The solid line is obtained with J = 142 K and a = 0.24. 

FIG. 3. A comparison between the experimental susceptibility of CuGeOs (filled squares, J2(| ) 
and results obtained from a 2D spin wave theory with J c = 115 K and = 0.1J C . 

FIG. 4. The two-magnon Raman continuum seen in experiment (filled squares, |32j, see also 
p9| ) at T = 20 K compared with results obtained from mean field theory (solid line). The 
theoretical results were obtained with J = 142 K and a = 0.24. 
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